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A systematic study of the root-mean-square distance between the constituents of weakly-bound 
nuclei consisting of two halo neutrons and a core is performed using a renormalized zero-range 
model. The radii are obtained from a universal scaling function that depends on the mass ratio 
of the neutron and the core, as well as on the nature of the subsystems, bound or virtual. Our 
calculations are qualitatively consistent with recent data for the neutron-neutron root-mean-square 
distance in the halo of ^^Li and ^''Be nuclei. 
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I. INTRODUCTION 

Light exotic-nuclei with one or two weakly bound neu- 
trons in their halo offer the opportunity to study large 
systems at small nuclear density. (A review on the theo- 
retical approaches and characteristics of halo nuclei can 
be found in ref. 0.) The constituents of the halo have 
a high probability to be found much beyond the interac- 
tion range. Then, the concept of a short-range interac- 
tion between the particles and its implications are useful 
in understanding the few-body physics of the halo. The 
quantum description of such large and weakly bound sys- 
tems are universal and can be defined by few physical 
scales despite the range and details of the pairwise in- 
teraction 131 . The particular halo- nuclei, with a neutron- 
neutron-core {n — n — c) structure, like ^He, ^^Li, ^^Be 
and ^°C, are examples of weakly-bound three-body sys- 
tems 0, where the above universal aspects can be ex- 
plored theoretically 

In weakly bound three-body systems, when the two- 
particle s-wave scattering lengths have large magnitudes, 
it is possible the occurrence of excited s-wave Efimov 
states 0. It was suggested in ref. that these states 
could be present in some halo nuclei. This possibility was 
further investigated and refined in ref. A parametric 
region was determined in which Efimov states can exist. 
Such region, for a bound three-body system, was defined 
in a plane given by the two possible (bound or virtual) 
subsystem energies. The promising candidate to exhibit 
an excited Efimov state was found to be [10. 

A few-body system interacting through a short raiige 
force can be parameterized by few physical scales |8|. 
For a zero-range force in three-space dimensions, it is ex- 
pected a new physical scale for each new particle added 



to the system, unless symmetry and/or angular momen- 
tum forbids the particles to be near each other. The 
three-body system in the state of zero total angular mo- 
mentum, has the bound or virtual subsystems energies 
and the ground state three-body energy as the dominat- 
ing physical scales. Any observable can be expressed in 
terms of the ratios between the physical scales, given by 
a scaling function, when the scattering length goes to 
infinity, or the interaction range tends to zero (scaling 
limit) In that sense the scaling function is an 

useful tool to study three-body observables and provides 
a zero order approximation to guide realistic calculations 
with short-range interactions. 

The scaling functions allow one to easily perform sys- 
tematic studies of various three-body halo-nuclei proper- 
ties with different types of two-bo dy s ubsystems as, for 
example, ^He, ^^Li, ^''Be and jlj. A classification 
scheme proposed for a bound three-body system j can 
be investigated in terms of a scaling function for the dif- 
ferent average distances between the constituents of a 
neutron-neutron-core halo nucleus. The classification of 
the three-body halo system depends on the kind of the 
two-body subsystems, bound or virtual. 

For the case of two identical particles in the three- 
body system, we have to consider four possibilities for the 
two-body subsystems (see figure 1): all unbound {Bor- 
romean); all bound {All-hound); one bound and two un- 
bound ( Tango configuration) ; and one unbound with 
two bound (we suggest a name Samba for this configura- 
tion). One natural example of halo- nuclei Samba system 
is the nucleus, which is composed by two-neutrons 
and a core. The neutron and forms the weakly 
bound state of ^^C 

In the present work, we study the root-mean-square 
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FIG. 1: Diagrammatical representation of the classification 
scheme for weakly-bound three-body systems. Two-body 
bound state are represented with solid line, and virtual state 
with dashed line. 



method to treat the Faddeev equations with two iden- 
tical particles, leading to the form factors from which we 
obtain the different mean-square radii. We also give a 
brief discussion of the scahng functions to describe the 
radii. In section III, we discuss the classification scheme. 
In section IV, we present our numerical results for the 
root-mean-square distances between the particles. Our 
conclusions are summarized in section V. 



distances between the constituents of a bound three-body 
system, where we have two identical particles n and a 
core named A. By n we mean neutron in halo nuclei, 
but we allow the pair neutron-neutron (nn) to be bound 
in order to cover all the configurations presented above. 
We represent the radii in a scaling plot in terms of a di- 
mensionless product, extending to halo nuclei a previous 
application that was done for molecules [ISi] . Using these 
scaling plots we can follow the behavior of the different 
radii when it happens a transition between one kind of 
system to another one. Starting from the Borromean case 
(all unbound), we can go to the Samba case by increasing 
the binding energy of the pair nA (keeping nn unbound); 
and to the Tango case by increasing the binding energy 
of the pair nn (keeping nA unbound). In particular, we 
calculate the mean square distance between the two neu- 
trons in ^He, ^^Li and ^^Be (which have been measured 
recently [T^ l from the known bound or virtual energies 
of the subsystems and the two-neutron separation energy. 

The calculation of the scaling function for the differ- 
ent radii of the three-body system is performed with a 
renormalization scheme applied to three-body equations 
with s— wave zero-range pairwise potential, which makes 
use of subtracted equations H, [13 • 

As our approach is restricted to s— wave two-body in- 
teractions, some limitations in the application of our 
analysis are expected. The argument has the origin in 
the fact that, in some cases, the two-body subsystem in- 
teraction in p— wave is considered to be important for the 
three-body binding. The examples are ^He known 
to be bound by the p— wave interaction in n—*Tie; and 
"'^^Li, where both p— and s— waves might be relevant 
to describe the ground state of the unbound subsystem 
^°Li (see discussion and references in 2^). However, as 
pointed out in ref.'2^ (in a discussion related to ^^Li), 
one should also noticed that even the three-body wave 
function with s— wave nn correlation produces a ground 
state of the halo nuclei with two or more shell-model con- 
figurations. The effect of the above s— wave restriction is 
also reduced due to the fact that we are considering the 
experimental energies in our approach, which implicitly 
are carrying the effect of higher partial waves in the in- 
teraction. Another aspect is that the radius are obtained 
from the tail of the wave-function, which is dominated 
by the s— wave. 

The paper is organized as follows. In section II, we 
present our formalism, which contains the subtracted 



II. FORMALISM 

In the next subsection, we consider the formalism given 
in refs. 0, for coupled channels to calculate the 
three-body wave functions and the possible different radii 
with a zero-range pairwise interaction. We solve the ho- 
mogeneous three-body Faddeev equations for a system 
with two identical particles a and a third one /3 in a 
renormalized subtracted form, which allows one to obtain 
the observables as a function of the two and three-body 
scales of the system. The corresponding masses of the 
particles a and /3 are nia and mp. 



A. Subtracted Faddeev Equations 

Next, we follow the model presented in refs. [isl llTj 
for the subtracted Faddeev equations, and consider units 
such that h = 1 and TOq = 1. For the subtraction energy 
that is required in the model, we choose /i^g^ . In this case, 
all the energies and momentum variables are rescaled to 
dimensionless variables considering this subtraction en- 
ergy. After partial wave projection, the s— wave coupled 
subtracted integral equations for the Faddeev spectator 
functions Xaa and Xa[3, are given by: 
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where we are defining the mass ratio and the constant 
mass factors by 



A = 



1 

2A 



(7) 



The plus and minus signs in eqs. (|3J| and Q refer to 
virtual and bound two-body subsystems, respectively. In 
the equations above, the dimensionless three-body energy 
£3 and the two-body energies (e^a and Eq/j), are related to 
the corresponding physical quantities by £3 = —E^/^'^^^, 

taa = -Eaal and tap = -Eap/fJ-f^y The thrce- 
body physical quantities can be written in terms of a 
scaling function, i.e., the dimensionless product of the 
observable and some power of the three-body binding 
energy E3, when the value of /i^g^ is determined from 
the known value of i?3 and consequently the three-body 
quantities are naturally a function of E3 and the ratios 
Eaa/Es and Eap/E^. Note that we are using the same 
symbol A for the mass ratio as well as for the core label, 
considering that, in a n—n— core nucleus we have the core 
consisting of A nucleons and A can also be identified with 
the mass ratio = Arun, rua = mn)- However, our 
expression can be generally applied for non-integer values 
of ^. 

For large scattering lengths, the details of the interac- 
tion are unimportant to describe few-body systems, as 
the short-range informations, beyond the two-body scat- 
tering lengths, are carried out by one typical three-body 
scale. Therefore, the low-energy observables present 
a scaling behavior quite universal with the three-body 
binding energy |^ |^ . In the limit of infinite scattering 
lengths or zero range interactions, the function which rep- 
resents a given correlation between two three-body ob- 
servables, written in terms of scaled variables, converges 
to a single curve, despite the existence of many other 
Efimov states. 

For a three-body system with binding energy E3, in 
the scaling limit |j, , one general three-body physical 
observable O, with dimension of energy to the power rj, 
at a particular energy E, can be written as a function J- 
of the ratios between the two and three-body energies, 
such that 



0{E,E3,E. 




(8) 



The two-body energies E'^^ (7 — a, /?), are negative 
quantities, corresponding to bound or virtual states. The 
nature of such two-body state, bound or virtual, is re- 
vealed in the momentum space, such that we have a 
bound state when \ Eap \ is positive and a virtual state 
when yJ\Eai3\ is negative. So, in equation ||SJ|, the signs 
-I- or — mean a bound or virtual two-body subsystem, 
respectively. The different radii of the bound aa/? sys- 
tem are functions defined from the eq. © with E = E^, 
which depend on the mass ratio, A, the ratios of the two 



and three-body energies and the kind of subsystem inter- 
actions (bound or virtual). 



B. Radii calculation 

The Faddeev components of the wave-function for the 
aa(3 system are written in terms of the spectator func- 
tions, obtained from the solution of eqs. iCJ and 101 m 
momentum space: 
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where is defined in eq. iQ. The Faddeev components 
are denoted by the sub-indices of the interacting pair. 
Representing the pair by aj with 7 = a or /5, one has 
that z is the relative momentum between a and 7 and 
y is the relative momentum of the third particle to the 
center-of-mass of the system aj. For the halo nuclei the 
notation is a = n and f3 is the core represented by A. 
The momentum component of the total wave-function. 



(10) 



is symmetrical by the exchange between the neutrons, 
n and n', while the antisymmetry is given by the sin- 
glet spin-component (not explicitly shown). The differ- 
ent mean-square radii are calculated from the derivative 
of the Fourier transform of the respective matter density 
in respect to the square of the momentum transfer. The 
relative mean-square distances between the halo neutrons 
and between the neutron and the core (7 — n. A) are ob- 
tained from the expression 



(11) 



where 



Fnjiq ) = I d^yd-'z^Ann' {y,z + ^] *a„„- 



y,^^-f ()i2) 



is the Fourier transform of the two-body densities, which 
is a function of the momentum transfer q (given in units 
of /i(3)). The relative momentum between n and 7 is 
z ± q/2 and the relative momentum of the third particle 
to the center-of-mass of wy is y. 

Analogous equations can be found for the mean square 
distances of the neutron and the core to the center-of- 
mass system ((r^)^), in terms of the Fourier transform of 
the one-body density. 
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C. Scaling functions for the radii 

The scaling functions for the mean-square separation 
distances between the particles in the three-body system 
can be written according to eq.||SJ) with E = E^. The 
scaling functions for the radii are written as: 

vW^=^r(±/^'±/^-^)- (14) 

To study the different types of three-body systems, the 
general scaling function for the radii given by eqs. l(T^ 
and H14|l will be studied in the configurations of figure 1 . 
However, one noticeable situation occurs when the Efi- 
mov limit is reached, for which the subsystems energies 
vanishes, and 

^{rl^)\E,\^nn, {A), y'(r^)|i?3|=7e7(A), (15) 

depends only on the mass ratio A (7 = n,A). Curi- 
ously, this configuration contains simultaneously all types 
shown in figure 1. 

III. CLASSIFICATION SCHEME: 
QUALITATIVE PROPERTIES 

A discussion of a classification scheme for a bound 
three-body system, with two identical particles, is given 
in ref . , according to the nature of the subsystem in- 
teractions, which can present a bound or virtual state, as 
depicted diagrammatically in figure 1. 

The different possibilities of three-body systems are 
reflected in the qualitative properties of the dynamics as 
given by the coupled equations 1^ and (O in terms of the 
strength of the attractive kernel of these equations. Let 
us describe all these possibilities. The Borromean case 
corresponds to positive signs in front of the square-root 
of the energies of the subsystems in both eqs. (|3J) and 
Q . A Tango three-body system has a negative sign only 
in front of ^/eaa in eq. JSJ , with positive sign in front of 
^eap ill eq. Q . For a Samba configuration of the three- 
body system, a minus sign appears multiplying JT^i in 
eq. while a plus sign multiplies ^taa in eq. (O. The 
All-Bound system has negative signs in eqs.© and 
As a consequence of the differences in eqs.® and the 
weakest attractive kernel in eqs.([Pl and is given by 
a Borromean three-body system, followed by the Tango, 
Samba and All-Bound systems. Therefore, once fixed a 
three-body binding energy, the system for which the ker- 
nel presents the weakest attraction will have the smallest 
configuration. So, the size of the corresponding system 
will increase in the following order: Borromean, Tango, 
Samba and All-Bound. The All-Bound configuration has 
the biggest size among all, for a given three-body binding 



energy. One of course has to remind that the sizes are 
expected to increase when the binding energy hits one 
scattering threshold. Therefore, for nonvanishing three- 
body binding energy this situation does not happen only 
in the Borromean case. In this sense, it is physically sen- 
sible that the Borromean case corresponds to the smallest 
configuration size. This observation will be explored in 
our numerical calculations. In this respect, we are show- 
ing the quantitative detailed implication to the different 
radii of weakly-bound three-body systems of the classifi- 
cation scheme proposed in ref. p^. 



IV. NUMERICAL RESULTS FOR THE RADII: 
SCALING PLOTS 

We present numerical results for the different possible 
radii for Borromean, Tango, Samba and All-Bound three- 
body configurations obtained from the wave-function, 
eq. l(Tn|) derived from the solution of eqs. ^ and In 
particular, we give results for the neutron-neutron (nn) 
average distance in the three-body halo nuclei, ^He, ^^Li, 
i^Be, 20c. 

It is interesting to present results in a general scaled 
form in terms of scaled two-body binding energies and 
mass ratio, as given by eqs. and 114|) . which turns to 
be useful for the prediction of the several radii in different 
situations for weakly bound molecules up to light halo- 
nuclei. Our calculations present results independent on 
the detailed form of the interactions in a weakly bound 
three-body system. It means that they apply very well 
to halo nuclei and weakly bound molecules. The present 
approach is valid as long as the interactions within the 
three-body system are of short range while the two-body 
energies are close to zero, i.e., the ratio between the in- 
teraction range and the modulus of the scattering lengths 
should be somewhat less than 1. These are indeed the 
cases we are considering. We present results only for 
the ground state of the three-body system. We have 
shown that the scaling function for the radii is indeed 
approached even for a calculation of the ground state in 
our subtraction scheme [15| , which will be enough for our 
purpose. 

In figure 2, we show the results of the scaling func- 
tion for the mean square distances, ^ (rn-y) and y^(r^) 

with 7 = n or A (see eq. ()15|l') which are shown as a 
function of A for Enj = 0. The comparison with ex- 
perimental results of \/{r^^ [l^ . for the ^He, I'^Be and 
^^Li is just for illustrative purpose, considering the hypo- 
thetical cases that ^He, ^'^Be and ^"Li would have virtual 
states close to zero energy. Such hypothesis is more re- 
alistic in case of ^"'Be, while for the ^He and ^^Li are 
not so. As shown by our results for ^^Li, given in Table 
I, the assumption of a virtual state with energy close to 
zero i\EnA\ ^0.05 MeV 20]) for lead to an average 
nn separation distance (in the halo of ^^Li) not compat- 
ible with the corresponding experimental result. Given 
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FIG. 2: Dimensionless products {rn^)\E3\ and (r^)|ii'3 
(7 = n, A) as a function of A for zero two-body energies. Re- 
sults for; (a) ^ {r?,„)\E3\ (upper solid line) and VvIa)T^ 
(upper dot-dashed line); (b) {r^)\E3\ (lower solid curve) 
and (r^) jiSg | (lower dot-dashed curve). The experimental 
data, obtained from are for -y/ (r^n)i-E3|- 



that it is well documented the difficulty in studying the 

we can consider ^ (r^n) ^ O'^^ of '^^^ 
puts to predict the virtual state of the n— ^Li system. In 
this case, we conclude that it cannot be smaller than ~ 
0.1 MeV. 

In this case {En^y — 0) the only physical scale is E3 
and the scaling plots will depend solely on A. The av- 
erage separation distance between the neutrons and the 
neutron-core tends naturally to a constant for large A, 
while it diverges for A tending towards zero. The reason 
for the unbound increase of the products {rnn)\^3\ ^^d 
(r^^)|i?3| for small A is due to the average momentum 
of the core which tends to zero (~ -^/^jiJal) extending 
the system to infinity. We have checked that the results 
shown in figure 2 tends to finite values after multiplica- 
tion by y/A (one has to remember that we are using units 
of m„ = 1). 

The average distance of the neutron to the center of 
mass system tends to become the relative distance to 
the core, when the core mass grows to infinity. This 
fact is clearly seen in figure 2 with the lower solid line 
approaching the upper dot-dashed one when A > 10. 
Also, one sees that the core distance to the center of 
mass vanishes with growing A as it should be. 

In the figures 3 to 6, we show results for VvmJP?3| 
and (r^^)|-E3|, when one of the subsystem energies is 
fixed in respect to the three-body binding energy, while 
the other one varies. We use values of A equal to 0.1, 1 
and 200. The subsystems can be bound or virtual and 
therefore all configurations are covered, i.e., we show re- 



TABLE I: Results of the neutron-neutron root-mean-square 
radii in halo nuclei. The cores (A) are given in the first col- 
umn, the absolute values of the three-body ground state en- 
ergies E3 are given in the second column. For —£3, which 
is equal to the two-neutron separation energy 5(271), we con- 
sider the center value of the corresponding energies given in 
ref. (3, except for Lithium. In case of Lithium we consider 
for —E3 the maximum (0.32 MeV) and the center value (0.29 
MeV) given in ref. [2^. In the third column we give our in- 
put values for —En a, considering several values, covering the 
values suggested in the literature (the references are given to- 
gether with the corresponding number) . For bound two-body 
subsystem nA, we have —EnA equal to the one-neutron sep- 
aration energy S{n). The virtual states are indicated by (v), 
and the nn virtual state energy is taken as -143 keV. The 
experimental values, in the last column, are obtained from 
Marques et al. (2000) 
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suits for Borromean, Tango, Samba and All-bound-type 
systems. Anticipating the presentation of these figures, 
in general we find that, the radii increase in the follow- 
ing order Borromean, Tango, Samba and All-bound for 
a given three-body energy and fixed A. In our analysis 
below, we fix either Enn/E^ = 0.1 or Eua/E^ = 0.1 
which can correspond to bound or virtual subsystem 
energies. In the next, we also consider the definitions 
EnA/E^ = Kf^j^ and Enn/E^ = K"^^ (for both bound or 
virtual-state energies), such that 



KnA = ±^EnAlE3 and Knn ^ ±VeZJe'3. (16) 

The +(— ) sign refers to bound (virtual) state. One should 
also note that, the dimensionless quantities Knn and KnA 
are directly related to poles in the imaginary axis of the 
respective two-body momentum space. 

In figure 3, we present calculations of the nA and 
nn root-mean-square radius as functions of KnA/\Knn\- 
Such average radius, multiplied by -y/j^Eaj, are scaled 
to dimensionless quantities. We consider A — 0.1 and 
Enn/Ej, = 0.1 fixed, corresponding to bound {Knn = 
VO.l) or virtual {Knn = — VO.l) subsystems. In this 
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FIG. 3: Dimensionless products \/ ('■^a)!-^^! (upper frame) 
and ^ (r^n)|-E3| (lower frame) for A = 0.1 and Enn/Ez — 
Knn = 0.1 as a function of KnA/\K„n\- In the upper frame, a 
bound nn pair is represented with dashed line; and a virtual 
nn pair with dot-dashed line. In the lower frame, a bound nn 
pair is represented with solid line; and a virtual nn pair with 
dotted line. The transition between the configurations occurs 
when KnA ~ (represented by the vertical line). 



case, the mass of the particle n is much heavier than 
the "core" , which does not happen in halo nuclei. How- 
ever, we consider this case for the sake of general in- 
terest. In the a;— axis, the positive (negative) values of 
KnA/\Knn\ correspond to bound (virtual) nA states. In 
the upper frame, the average nA radius is shown for a 
bound nn pair (dashed line) and for a virtual nn pair 
(dot-dashed line). So, the dashed line shows that the 
value of \/ (r'nA) I -^3 1 increases with the transition from 
Tango to All-bound configuration. In the other possibil- 
ity, represented by the dot-dashed line {nn in a virtual 
state), the value of \/ (^^a) 1-^3 1 increases from the most 
compact Borromean configuration {KnA negative) to the 
Samba-type configuration {KnA positive). 

In the lo wer frame of figure 3, we also show that 
■y/ (?'^i,i)|-E'3| increases with the transition from Tango to 
All-bound and Borromean to Samba configurations. We 
observe, in this case, a small sensitivity on Knn, when 
going from — \/0.1 (dotted line) to VO.l (solid line), with 
\/Jr^m)\^3\ having practically the same value for the All- 
bound and Samba configurations; and also for the Tango 
and Borromean configurations. As we increase A, this 
sensitivity increases, as one can see in the next figures 4 
and 5. We interpret this as the dominant role played by 
a light particle in the long-range interaction between two 
heavy particles, n, as already shown in an adiabatical 
approximation of the three-body system . 

In figure 4, we show the radii for the ground state of 
the three-body system with the mass ratio A = 1. As 
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FIG. 4: Dimensionless products ^ {^\ai\^^\ (upper frame) 

and -y/ {r'^n)\Ei,\ (lower frame) for A — 1 and — 0.1 as 
a function of KnA/\Knri\. The convention of the lines is the 
same as given in figure 3. 
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FIG. 5: Dimensionless products \/ ('"^a) 1-^3 1 (upper frame) 
and ^y {r'^„)\E3 \ (lower frame) for A — 200 and Enn/Ez = 0.1 
as a function of KnA/\Knn\- The convention of the lines is 
the same as given in figures 3 and 4. 



in figure 3, we fixed Knn — ±\/0.1, corresponding to 
bound (+) or virtual (— ) subsystems. The same behav- 
ior found in figure 3 for VvIVT^^ i^ found in the upper 
frame of figure 4 for A = 1, i.e., the configurations for 
which the nn pair is virtual (dot-dashed line) are smaller 
than the ones that have the nn pair bound (dashed line) . 
Besides that, the configuration increases in size when 
the system changes from a Borromean to a Samba type 
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FIG. 6: Dimensionless products vvLi)T^^ (upper frame) 
and ^y{rl„)\E3\ (lower frame) for A = 200 and E„a/E3 =0.1 
as a function of Knn/\K„A\- In the upper we represent the 
case of a bound nA pair with dashed line; and the case of a 
virtual nA pair with dot-dashed line. In the lower frame we 
represent the bound nA pair with solid line and the virtual nA 
pair with dotted line. The lines of this figure correspond to 
the vertical transitions represented in figure 5 in both extreme 
sides where l/i'nAl ~ \Knn\- 



and when it changes from a Tango to an All-bound type. 
For A — 1, the mean square distance between the two- 
neutrons (lower frame of figure 4) exhibits the same qual- 
itative behavior as found for \/ ('^^^1)1-^31 when the con- 
figuration type is modified for a fixed three-body energy. 
These conclusions are still valid for a heavy core particle 
with A = 200, as one can verify in figure 5. It is worth- 
while to mention that \/ 1-^3 1 attains its asymptotic 
value fast with the increase of A consistently with the 
calculations presented in figure 2 for = EnA = 0- 
(Compare the results for A = 1 and A = 200 in the lower 
frames of figures 4 and 5). 

In correspondence with figure 5, for the same mass ra- 
tio A = 200, in the last set of systematic calculations 
we consider in figure 6 the energy of the subsystem nA 
fixed in relation to E^, with KnA = ±V0.1, correspond- 
ing to bound or virtual (— ) nA subsystem. In this 
case, the ratio Knn/\KnA\ is changed from —1 to +1, 
such that the transitions of configurations from the left 
side to the right side of this figure correspond to the 
vertical transitions in the extreme side of figure 5 (when 
= \Knn\ = VO.l). So, comparing the upper frames 
of figures 5 and 6, we observe that the vertical transition 
from All-bound to Samba configuration in figure 5 corre- 
sponds to the dashed line of figure 6; and, the vertical 
transition from Tango to Borromean configuration in fig- 
ure 5 corresponds to the dot-dashed line of figure 6. In 
the lower frames of both figures, similarly we observe 



that the vertical transitions of figure 5 correspond to the 
lines represented in figure 6: the transition from Samba to 
All-bound configuration is given by the solid line; and the 
transition from Borromean to Tango configuration given 
by the dotted line. In general, one can observe that the 
three-body bound state has a smaller size when the nA 
pair forms a virtual state (see in each frame of figure 6, 
where the upper line is for bound and the lower line is 
for virtual nA pair). 

The calculations of the average distances of the neu- 
trons in the halo of *^He, "Li, "Be and ^^C are shown 
in Table I and compared with the available experimental 
data. For the input, we have considered i?„„ — —0.143 
MeV and the center of the available experimental values 
of E^, and EnA- For the cases that we have unbound vir- 
tual nA systems, the EnA input values are taken from 
several recent theoretical and experimental analysis; we 
prefer to keep at least two possible values in order to 
verify the consistency of the model with the experimen- 
tal data. 

Within the possible limitations of our approach, by 
comparing our results for the nn mean-square radius with 
the experimental ones, which are known in the cases of 
®He, ^^Li and ^^Be, as given in Table 1, one can also 
predict the corresponding virtual energies of the nA sub- 
system. 

In the particular case of ^He, the comparison between 
the experimental nn mean-square radius with our result 
is pointing out to a virtual energy close to zero for n— '*He, 
which is not supported by the quite large values given in 
the literature [ll, . However, as discussed in |3 , 
the interaction for n— ^He is attractive in p— wave and 
repulsive in s— wave producing a large value for the en- 
ergy of the virtual state, such that a deviation of our 
calculation from the experimental result is expected (in 
our model the s— wave poles should be near the scatter- 
ing thresholds). For the binding of ^^Li, it is also known 
that both p— and s— waves are important in the subsys- 
tem ^Li. This fact can also explain some deviations 
of our results when compared with experimental ones. 
However, in this last case, our approach can be more 
reliable based on the fact that: (i) the s— wave n—^Li 
interaction is attractive and it has a virtual state near 
the scattering region; (ii) we are considering the experi- 
mental energies for the inputs, such that we are partially 
taking into account the effect of higher partial waves in 
the interactions; (Hi) the radius is obtained from the tail 
of the wave-function, which is dominated by the s— wave; 
(iv) as noted in ref. (2ft] . even the three-body wave func- 
tion with s— wave nn correlation produces a ground state 
of the halo nuclei with two or more shell-model configu- 
rations. 

One should also expect that other effects like the finite 
size of the core and Pauli principle, missing in our model, 
could affect the average relative distances, if the three- 
body wave function would overlap appreciably with the 
core. At least for ^^Li this is not expected due to the 
small binding. It is of notice that indeed the halo neu- 
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trons have a large probability to be in a region in which 
the wave function is an eigenstate of the free Hamilto- 
nian, and thus dominated by few asymptotic scales. 

V. CONCLUSIONS 

The mean-square radii of a light halo-nuclei modelled 
as a three-body system (two neutrons n and a core A) are 
calculated using a renormalized three-body model with a 
pairwise Dirac-5 interaction, which works with a minimal 
number of physical inputs directly related to observables. 
These physical scales are the two-neutron separation en- 
ergy S{2n) — and the nn and n— core s— wave scat- 
tering lengths. 

The existent data for ^^Li and ^''Be compare qualita- 
tively well with our theoretical results, which means that 
the neutrons of the halo have a large probability to be 
found outside the interaction range. Therefore the low- 
energy properties of these halo neutrons are, to a large 
extend, model independent as long as few physical input 
scales are fixed. The model provides a good insight into 
the three-body structure of halo nuclei, even considering 
some of its limitations. We pointed out that the model 
is restricted to s— wave subsystems, with small energies 
for the bound or virtual states. So, the model is not suit- 
able for the ®He, since the s— wave virtual state energy of 
^He is quite large (-4 MeV). There is no n— core p— wave 
interaction, although some of its physics is effectively in- 
cluded through the value of the two-neutron separation 
energy, which is an input for our radii calculations. Also 
the finite size of the core and consequently the antisym- 
metrization of the total nuclear wave function, are both 
missing in our model. However, as the three-body halo 
nuclei tend to be more and more weakly bound with sub- 
systems that have bound or virtual states near the scat- 
tering threshold, our approach becomes adequated and 
the above limitations are softened. The results indicate 
that the model is reasonable for ^^Li and ^''Be. 

As an example of application to other halo-nuclei sys- 
tem, considering the available energies, we have also esti- 



mated the nn root-mean-square radius for the n — n—^^C 
system. 

We have also studied in detail the consequences of the 
classification scheme proposed in Ref. [l^ for weakly 
bound three-body systems. This study was performed 
analyzing the dimensionless products ^/Jt^^^a^WsI ^^'^ 
V (^nn)l-^3l ™ terms of scaling functions depending on 
the dimensionless product of the scattering lengths and 
the square-root of the neutron-neutron separation energy. 
In the cases we have addressed, there are four different 
types of a three-body system when we allow the neutron- 
neutron pair to be bound: Borromean (only virtual two- 
body subsystems). Tango (nn bound and nA virtual). 
Samba (nn virtual and nA bound) and All-Bound (only 
bound two-body subsystems). The name Samba was in- 
troduced to refer to a system quite stable because it has 
two bound two-body subsystem than the Tango type, so 
it can "shake" a little bit more and continue to be bound. 

The qualitative properties of the different possibilities 
of three-body systems are easily understood in terms of 
the effective attraction in the model: when a pair has 
a virtual state the effective interaction is weaker than 
when the pair is bound. Thus, a three-body system has 
to shrink to keep the binding energy unchanged if a pair 
which is bound turns to be virtual. We have illustrated 
through several examples this property, which show that 
dimensionless sizes VvIa)T^^ ^^'^ VvmiTP^ increase 
from Borromean, Tango, Samba and to All-Bound con- 
figurations. Of course the size is expected to increase 
beyond limits when a nonvanishing three-body energy 
hits a scattering threshold, with the Borromean config- 
uration being the only exception. In spite of that, we 
conclude that even far from the threshold situation, the 
configuration sizes increase as we pointed out. 
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